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1. Understanding the modification of hadron properties in a hot and dense medium is one of the most important 
issues, since it is deeply related to the restoration of chiral symmetry of quantum chromodynamics (QCD) 0, Q. In 
particular, it is of great importance to understand the change of hadron masses in medium, since much of them is 
generated by the spontaneous breakdown of chiral symmetry (SB%S) of QCD. The quark condensate is regarded as an 
order parameter of SB^S and is known to decrease as temperature and density increase. Eventually chiral symmetry 
will be restored at certain critical temperature and density. Thus, it implies that the masses of hadrons are expected 
to be modified in medium. 

The mass dropping of the vector mesons in dense matter has been extensively studied theoretically as well as 
experimentally (see a recent review Q and references therein). In the early 1990s, Brown and Rho investigated 
the mass dropping of hadrons in dense medium, based on the effective chiral Lagrangian with scale invariance, and 
suggested the in-medium scaling law, also known as the BR scaling According to the proposed scaling law, the 
masses of the vector mesons are expected to drop approximately 20 % at the normal nuclear matter density. Hatsuda 
and Lee drew a similar conclusion within the QCD sum rules [f|. A most recent reanalysis of the QCD sum rules 
suggests about 15 % dropping of the p meson mass at the normal nuclear matter density [f|. The medium modification 
of the p and co mesons in cold nuclear matter was experimentally reported only very recently. The KEK-PS E325 
experiment measured the invariant mass spectra of p, u> — > e + e~ decay modes in proton- induced nuclear reactions Q 
and announced that the masses of the p and ui mesons decrease by about 9 % at the normal nuclear matter density. 
On the contrary, the CLAS experiment found almost no evidence of the mass dropping of the p meson [8, 9]. Instead, 
the CLAS experiment observed the broadening of the p meson decay width that is in line with the predictions by 
hadronic models [Io|, [H| • 

The nucleon in medium also undergoes modification and eventually its structure gets altered. A typical example is 
the well known nuclear EMC effect [12(. A recent experiment of the ratio of the electric and magnetic form factors 
of the proton in the 4 He(e, e'p) 3 H reaction indicates that the electromagnetic (EM) form factors are changed in 
nuclei 1131 . The modification of the nucleon mass in nuclear matter has been studied within various nuclear models 



such as the Walecka model (see a review [14|, and references therein), the quark- meson coupling (QMC) model (see a 
review [lj|, and references therein), the soliton-bag model [lij], the Skyrmc model [13, (3 and so on. However, there 
is no consensus about how much its mass decreases so far. Depending on models, the effective nucleon mass (M^) at 
the normal nuclear matter density is known to be dropped in the range of 500 MeV < < 930 MeV. On the other 
hand, Celenza et al. fl9j proposed that the size of the nucleon increased in nuclear medium in order to explain the 
EMC effect. 

In this Letter, we investigate the bulk properties of symmetric nuclear matter within the framework of a chiral 
soliton model based on the effective Lagrangian with the n, p, and u meson degrees of freedom. This chiral solitonic 
approach has a certain virtue that the mesonic sector governs the baryonic sector, so that the medium modification 
of the mesons and nucleon can be treated on the same footing. The chiral soliton model with vector mesons in 
free space [20l - [2"i ] has been developed in 1980s such that the p and u> mesons stabilize the skyrmion in place of the 
original stabilizing Skyrme term [25T - [27| . The Skyrme term is known to be related to the massive vector mesons [28| . 
The chiral soliton model with vector mesons was shown to describe well the static properties of the nucleon in free 
space (see the review [24| and references therein). In the meanwhile, the original Skyrme model was extended to 
finite density [l7j . the pion properties being modified in a dense medium based on the experimental data of pion- 
nucleus scattering. In Ref. [1 81 ] , the model was further elaborated by revising the stabilizing Skyrme term in nuclear 
matter. The medium-modified Skyrmc model improved the selfconsistcncy of the approach and could describe the 
basic properties of nuclear matter. In the present work, we propose an in-medium chiral soliton model with vector 
mesons. The model has an advantage that it describes the meson and nucleon properties in free space and in nuclear 
matter, and the bulk properties of nuclear matter on an equal footing. In particular, the model shows one possible 
way as to how recent issues of the mass dropping of the vector mesons can be related to the properties of nuclear 
matter. 

2. Our starting point is the effective chiral Lagrangian with the n, p, and u> meson degrees of freedom [22l. |23T| . where 
the nucleon arises as a topological solution. We extend the model by introducing the medium functionals that encode 
the influence of the surrounding nuclear environment to the nucleon structure and the corresponding properties in 
line with the previous works [13, [H| , where the "outer shell" [TtJ and "inner core" [lH modifications of the nucleon 
in nuclear matter have been considered. The corresponding in-medium modified Lagrangian has the following foruQ 

C* = + Cy + C^y + C\y Z , (1) 



1 Hereafter, the asterisk indicates explicit medium modifications. 
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where the kinetic and mass terms are written as 

tl #2 



£* = 



f -f Tr (doUdotf) - a p (p) ^Tr {8 l Ud l U^) + a s (p) Tr (U + rf-2), 



c *v = -^-^^{W-d v v;-i[v;,v:]} , (2) 

where stands for the generic coupling constant for a vector meson. The vector fields arc defined as 

Vjj>> = 9 l T a pl, V<?> = f w„, (3) 
and the interaction part has the following form 

Cl v = 2f 2 Tr (j p V<f> Vj?>) 2 . (4) 



M ' fJ. 

The vector current with the it field is defined as 

~ 2i det(l + E/) ■ [ ' 

The last term in Eq.([T]) denotes the Wess-Zumino term 

£w* = \g>^ (6) 

which is related to the topological current 

B " = ^2^{(^dM)(U^d a U)(U^d,U)}. (7) 

The temporal component of Eq. ([7]) is identified as the density of the single baryon located in free space or in nuclear 
matter, satisfying J d 3 rB° = 1. 

The pion fields are parametrized as a form of the SU(2) matrix U = exp(zx • 7r// w ), while p° and w M represent 
the isovector-vector and the isoscalar-vector fields. We have introduced the additional medium functionals ( p (p) and 
Cui(p) such that the effective in- medium couplings are changed in nuclear matter as follows: g* = gCp(p) 1 ^ 2 an d 
9w = 9^{p) 1 ^ 2 - The input parameters /„■, to x and g stand for the pion decay constant, the pion mass and the ptttt 
coupling constant g = g plXT: in free space, respectively. The medium functionals a s (p) and a p {p) in Eq. {1]) are related 
to the phenomenology of low-energy pion- nucleus scattering [UHil- The ( p (p) and Cj(p) are considered in the kinetic 
term of the vector mesons and interaction parts phcnomenologically. Later, we will show that they will influence the 
inner core of the nucleoli. In fact, it is similar to the modification of the skyrme parameter e — > e* = eC, 1 ^ 2 as 
done in Ref. [l8[ . The modified Skyrme term presented in Ref. [l8| stabilizes not only the solitons of the model but 
also homogeneous nuclear matter. Similarly, the modified p and ui mesons stabilize the solitons and reproduce the 
properties of symmetric nuclear matter. We want to mention that the masses of the p and ui mesons are degenerate 
Skyri 

However, we consider two different models in the present work. First, we keep the masses of the p and lo mesons 
degenerate, assuming that the KSRF relation holds also in nuclear matter 



in the original Skyrme model with the vector mesons [22j, l23l ] and they are related to the input parameters by the 
KSRF relation j30[: ra 2 p = = 2f 2 g 2 in free space, i.e. in the case of a p (0) = 1, a s (0) = 1 and Cp(0) = Cj(0) = 1- 



2/^5 C = m v = m* = m* , C P = Cu 1 = C, (8) 

but the coupling constant in the Lagrangian will undergo a change in nuclear medium as g — > g v = gC, 1 ^ 2 . We call 
it Model I. On the other hand, we will also assume that the KSRF relation in nuclear matter is valid only for the p 
meson with the to meson intact, i.e. 

2/^ 2 C = m* V ml 2 = ml Q p = C, Cu, = 1. (9) 

This will be called Model II. In doing so, we can understand how the p and w mesons play roles differently in describing 
the properties of nuclear matter. The decay constant and mass of the pion in free space are fixed by the experimental 
data, i.e. f„ = 93 MeV and m ff = 135 MeV (the neutral pion mass). The coupling constant is chosen as g = 5.85, 
which is close to the empirical value .g p7r7r = 6.11 that is determined from the width of the 7 — >• 2-7T decay. The vector 
meson mass in free space is taken to be my = m p = m w = 770 MeV. 
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The medium functionals are given as [13, [l8, 29] 



1 



4irc Q p/r/ 
1 + g' ATrcop/r)' 



1 



cxp 



1 



7donP 



(10) 



where ?/ = 1 + 171^/171^ — 1.14 is a kincmatical factor. The empirical parameters bo = — 0.024 m~ x and cq = 0.09 m~ 3 
are in general consistent with the analysis of pionic atoms and the data of low-energy pion-nucleus scattering, and 
.9o = 0-7 is an empirical value of the correlation parameter [29j . 

Most solitonic models yield larger values of the nucleon mass and the present model gives tojv = 1526 MeV in free 
space whereas the mass splittings are correctly reproduced. Thus, instead of the absolute value of the nucleon 
mass, we consider the following parameter directly related to the coefficient of the volume term in the Weizsacker 
mass formula for nuclei as done in Ref. [la ]. One can minimize the mass functional derived from the Lagrangian 
and obtain the effective nucleon mass m* N in nuclear matter for a given set of parameters 7 num an d 7dcn- Then, the 



binding energy per nucleon can be defined as 
m* N - m N =e = -a v (p) H , 



(11) 



where ay(p) is the density-dependent coefficient for the volume term and the dots represent the surface, Coulomb, 
symmetry energy and pairing terms. Since we consider an infinite isospin symmetric matter, we concentrate only on 
the volume term and ignore all other terms in the nuclear mass formula. The density dependence of ay(p) is obvious, 
because it describes the binding energy per nucleon in symmetric nuclear matter. We rewrite Eq. (|ll[) as 



m N 

m N 



m N -a v {po) rn 



(exp) 
N 



,( gx p) 



(Po) 



m N 



,( cx p) 

N 



(12) 



so that the relation will be satisfied at the saturation point of nuclear matter, p = p$. The parameters 7 n 



(cxp) 
N 



939 MeV and 



and 7den are fitted to reproduce the correct saturation point. Note that experimentally m 

a,y Xp \po) = 16 MeV. However, almost all chiral soliton models produce the larger values of the nucleon mass. In the 
Skyrmc model with vector mesons, the nucleon mass typically turns out to be around 1500 MeV, while the A-nucleon 
mass difference in free space AM = Ma — Mjv is obtained to be similar to the experimental data. Since the nucleon 
mass is overestimated by about 40%, it is better to use av(po) = 26 MeV phenomenologically, which is about 40% 

larger than the typical value dy Xp \po) = 16 MeV. For example, in Ref. [lH, the experimental value of the nucleon 
mass, i.e. = 939 MeV was shown to reproduce ay = 16 MeV correctly. 

As mentioned already, we arc interested in symmetric nuclear matter. Since the nuclear density p is constant, one 
can use the spherically symmetric approximation for the meson profiles, i.e. the hedgehog Ansatzc 



U 



exp {it ■ nF(r)} , h 



Pi 



G(r) 
rL* 



w(r)<5 



(13) 



where the scale factor L* = g* is required to have the same boundary conditions for the medium-modified profile 
function G as in free space, i.e. G(0) = —2 (see Fig. [1]). 

Using the hedgehog Ansatze, one can derive the static energy functional from the Lagrangian, which is identified 
as the classical skyrmion mass 

M* H = E*[F(r) 7 G(r),w(r)} 

= - J £*d 3 r 

r°° ( f 2 

= 4tt / dr I a p ^ (r 2 F' 2 + 2 sin 2 F) + 2f 2 (G + 1 - cos F) 2 + a s r 2 f 2 ml (1 - cos F) 
G' 2 + ^t^ ] r 2 ^ 5 *V + \ 2 ) + g§ coF' sin 2 F \ , 



2r 2 



(14) 



where /' = df /dr, generically. Minimizing the classical soliton mass can be achieved by solving the equations of 
motion, which are given as the coupled nonlinear differential equations 



F" = --F' 



1 



G" 



r 

*2 f 2 



[4(G + 1) sin F — a.p sin 2F] + — ml sinF 



y. v r 2 1 1 u 1 ' a v kn 2 a p f 2 r 2 



— 7? sin^ F . 



2 



G + 2 cos ■ 



F 



G(G + l)(G + 2) 



J + 2flg*Juo 



3g: 

An 2 r 2 



F' sin 2 F 



(15) 
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with the boundary conditions 

F(0) = 7T, F(oo) = 0, G(0) = -2, G(oo)=0, o/(0) = 0, w(oo) = 0. (16) 
The collective quantization brings out the following relations 

U(r,t) = A(t)U(r)A+(t), Ui (f,t) = ^ (Kxn) t , 
T-p (r,t) = ^A(t)T-[K^(r)+n-(Kxn)Ur)]A + (t), 

r- Pi (r,t) = A{t)r- Pi {r)A + {t), (17) 

where 2K denotes the angular velocity of the soliton with the relation ir-K = A + A. This leads to the time-dependent 
collective Hamiltonian 

H*(t) = M* H {F,p,uj) - A*(F,p,Lu, ( p^ 1 ^ 2 )TiiAA + ), (18) 
where A* denotes the in-medium modified moment of inertia of the rotating skyrmion 

A* = g^F'si n 2 F + H/> 2 ^in 2 F + 8 S i n 4 ^-8e lS in 2 ^ + 2^ + (a+e 2 ) 2 ) 



S + 66 - 2^ - 6 + 1) + 2(G 2 + 2G + 2)£ + r 2 (3^f + g + 2&. 



p 



3.9, 

-\^ ,2 + ^)-\ f ^- (19) 

In the large jV c expansion, one extremizes the moment of inertia and gets the coupled nonlinear differential equations 
for the next-order profile functions £ij 6j 4> m the presence of the leading-order profile functions F, G and ui 

s = 2/^« 1 -i + cs F )-gi + G2 «-- 1 't 2(G+1)6 . 
« - W & + 1 - «^ - # + G ' fe -"' +2 r f +3G+3)& , 

,„ ,,2,2, 20 

* = 2/, 9 „ 4, — 2 + (20) 

with the boundary conditions 

0(0) = 0(oo) = 0, ei(0) - = 0, £a(0) = 6(00) = 0. (21) 

The boundary conditions for £1 and £2 satisfy the relation 2£i(0) + 6(0) = 2 that remains unchanged in nuclear 
matter. 

Finally, the effective masses of the nucleon and the A isobar are expressed in terms of the in-medium hedgehog 
mass Mfl and the in-medium modified moment of inertia A* 

™* N = M* H + ^, m * A = M* H + ^. (22) 

4. Figure [1] draws the solutions of the equations of motion (fT5|) for Model I, i.e. the profile functions for the 7r, 
p, uj background fields. The solid curves depict those in free space, whereas the dashed ones represent the modified 
profile functions at the normal nuclear matter density po = 0.15 fm~ 3 . One can clearly see that the profile functions 
F(r) and G{r) are enhanced everywhere. On the contrary, the uj one is suppressed up to around r ~ 0.6 fm, it starts 
to be strengthened. This is partly due to the fact that the profile functions F(r) and G(r) are constrained to be ir 
and 2 at r = by the Dirichlct boundary conditions, whereas uj(r) has the Cauchy boundary conditions. We also 
found that all the profile functions showed similar behaviors in the case of Model II. In Fig. [21 the profile functions 
£1, £2 and 4> f° r Model I are illustrated. These three profile functions also show enhancement in larger distances but 
are suppressed in smaller distances. The similar results are obtained also in the case of Model II. 

Figure[3]shows the results of the coefficient of the volume term ay as a function of nuclear matter density normalized 
by the normal nuclear matter density, i.e. p/po- The solid curve represents the result of Model I. In this case, ay 
increases as the density increases and then is already saturated to be approximately ay — 26 MeV at about p ~ 0.6 po 
and remains almost constant up to the normal nuclear matter density. Then, it decreases very slowly. It is very similar 
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FIG. 1. (Color online) Meson profile functions in free space (solid curves) and and at the normal nuclear matter density 
po = 0.15 fm -3 (dashed ones). The results are those from the Model I. 
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FIG. 2. (Color online) The profile functions £i, £2, and 10 <^>. The solid curves depict them in free space, whereas dashed ones 
represent them at the normal nuclear matter density po = 0.15 fm -3 . The results are those from the Model I. 



to those obtained within the in-medium chiral soliton model, where the Skyrme term was used instead of those with 
the vector mesons [3l|. This is due to the fact that the overall changes of the p and ui mesons in nuclear matter 
behave in almost the same way as the change of the Skyrme term as presented in Ref. (3~l| . On the other hand, we 
find that ay has a clear minimum for the binding energy at the normal nuclear matter density in the case of Model 
II, as shown in the dashed curve in Fig. [3] In Model II, the u meson remains intact, whereas the p meson undergoes 
a change in nuclear matter. It implies that a inner core part of the nucleon is less modified than its outer shell part. 
Thus, the result of Model II, which is more favorable (see the results in Table H]) than that of Model I in the present 
approach, can be interpreted as follows: as the density increases, the wavefunctions of nucleons in medium are forced 
to be overlapped but the effects of the Pauli blocking resist this overlapping. It means that while the outer shell part 
of a nucleon can be easily influenced, the inner core part is less affected by its neighborhoods. 

One can make a Taylor expansion of the coefficient of the volume term near the normal nuclear matter density with 
respect to p 

S a + I ^p- K )' + ^p- K f + .--, (23) 

where ao, Kq, and Qo denote the value of ay at p = po, the incompressibility of symmetric nuclear matter, and 
the coefficient of the third derivative, respectively. The coefficient of the first derivative, which is defined as the 
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FIG. 3. (Color online) The coefficient of the volume term in the binding energy formula (| 1 1 1) av as a function of nuclear matter 
density normalized by the normal nuclear matter density, i.e. p/po- The solid curve draws the result of Model I while the 
dashed one represents that of Model II. 

pressure, vanishes at the saturation point. In Table U we list the results of these three coefficients in comparison 
with other models. The parameters 7 num and 7d cn are fitted to reproduce ay ~ 26 MeV in both Model I and II. 
The incompressibility of symmetric nuclear matter Kq is an essential quantity, because it shapes the dependence of 





/num. 


/den 


Po [fm- 3 ] 


a [MeV] K [MeV] Q 


[MeV] 


m*/m p 


Model I 


2.390 


1.172 


0.150 


26.75 


110 


1591 


0.70 


Model II 


1.970 


0.841 


0.150 


26.38 


233 


1213 


0.72 


Hybrid [32J 






0.148 


16.24 


230 


-72 




FSU [33J 






0.148 


16.30 


230 


-523 




NL3 [34J 






0.148 


16.24 


272 


204 





TABLE I. The bulk properties of nuclear matter: po denotes the saturation point of nuclear matter density that corresponds to 
the normal nuclear matter density; ao represents the value of ay at the saturation point; Kq designates the compressibility of 
symmetric nuclear matter; Qo is the coefficient of the third derivative at po; m* P lm p is the ratio of the mass of the p meson at 
the normal nuclear matter density to that of the vacuum. The first and second rows list the corresponding results from Model 
I and Model II, respectively, while the other three rows are those from Refs. [32T - [34| | . 

the nuclear binding energy on the density near the saturation point. It also means that it shows how stiff or soft 
symmetric nuclear matter is. In the fifth column of TableU the results of the present models and other three different 
models are compared. Model I gives 110 MeV while Model II yields 233 MeV, which implies that the nuclear binding 
energy turns out to be stiffcr in the case of Model II, that is, it increases more rapidly with density than that of Model 
I does. The result of Model II is comparable with the predictions of the other models. Compared with the empirical 
value of the incompressibility extracted from the experimental data of the giant monopole resonance Kq = 235 ± 14 
MeV [HI , the present result of Model II is in good agreement with the data. The sixth column of Table Q] lists the 
predictions of Qo in comparison with those of the other models. The present results for Qo are quite larger than those 
of Refs. (32h34| . This is due to the fact that the medium functionals are determined from pion-nucleus scattering 
phenomenology at subnormal nuclear matter densities. If p is larger than p$, the effective parameters of the medium 
functionals and, in particular, the correlation parameter g' may become density-dependent. 

As mentioned earlier, an advantage of the present approach is that we can treat the mesonic and baryonic sectors on 
the same footing. Since we fix the relevant parameters to reproduce the properties of the nucleon and nuclear matter, 
we are able to find the p meson mass at a given density. In this way, we determine how the p meson mass is dropping as 
the density increases. Figure [4] depicts the mass dropping of the p meson as a function of the normalized nuclear matter 
density p/ po- In Model I, the masses of the p and u> mesons are degenerate and hence their changes in nuclear matter 
are identical. As shown in Fig. 21 they decrease as the density increases. However, their dependence on the nuclear 
matter density does not seem to be linear, i.e. it is not satisfied with the linear relation m p (p) = m p (0)(l — Cp/po) 
with constant C but it seems to behave as a quadratic function of the density. At the normal nuclear matter density 
po, the masses of the p and ui mesons are dropped by about 30% in the case of Model I. Model II brings about the 
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FIG. 4. The ratio of the mass of the p meson in nuclear medium to its value in free space as a function of normalized nuclear 
matter density p/po- The solid curve represents the result of Model I, while the dashed one draws that of Model II. 



p meson mass dropping by about 27%, while that of the u> meson is kept to be fixed. These results show a similar 
tendency to those obtained within the QCD sum rules @ and the generalized NJL model (36[, in which the mass of 
the p meson is approximately dropped by about 20%. 

As discussed previously, chiral soliton models yield in general the overestimated nucleon mass, so that we fix all the 
parameters to produce ay = 26 MeV, not the experimental value a^ p = 16 MeV. If the same argument applies to the 
mass dropping of the p meson, the 27 — 30 % mass dropping is more or less comparable with those of the QCD sum 
rules and NJL model, though it turns out to be larger than the recent measurement of the KEK-PS E325 experiment. 
However, we want to emphasize that in the present work we want to present a possible alternative way of explaining 
the changes of the nucleon and the bulk properties of nuclear matter from a different point of view. 



5. In the present work, we aimed at investigating the medium modification of the vector mesons and nucleon on 
an equal footing, and at studying the bulk properties of symmetric nuclear matter, based on the in-medium modified 
■n-p-u soliton model. We employed the two different models, i.e. Model I in which we keep the p and uj mesons 
degenerate, and Model II, where we let the p meson undergo a change while the u is kept to be intact. The pionic 
sector was determined by using pion-nucleus scattering data, while the vector mesons were assumed to be changed in a 
similar manner to the Skyrme term in Ref. [l8| . We first examined how the profile functions were changed in nuclear 
matter. Then, we determined the coefficient of the volume term in the nuclear mass formula. Since the nucleon 
mass from this model is overestimated, we also considered reasonably the larger value of the volume-term coefficient. 
We found that Model I produced the coefficient of the volume term that became almost saturated around ay — 26 
MeV. On the other hand, Model II exhibited a clear minimum of the binding energy per nucleon at around 26 MeV, 
from which the coefficient of the volume term started to decrease. The results from Model II were more favorable 
because of the following reason: as the density increases, the outer pionic shell of a nucleon is more influenced by its 
environment, while the change of its inner core is less recognized. We also calculated the incompressibility of nuclear 
matter, which is one of the most important bulk properties of nuclear matter. We obtained Ko = 233 MeV from 
Model II, which is in agreement with rclativistic mean-field models. Finally, we determined the mass dropping of the 
p meson that arose from the fixing of the vector-meson parameters. It was found that the p meson mass was dropped 
by about 27 - 30 %. 

The present work sheds light on the bulk properties of nuclear matter from a different point of view. While the 
soliton model with vector mesons provide a simple framework in investigating the properties of nuclear matter, it 
still brings about certain physical implications: in this framework, the modification of the vector mesons shows a 
certain correlation with the changes of the nucleon and nuclear matter properties. Moreover, the present approach 
can be easily extended to asymmetric nuclear matter or neutron-rich nuclear matter, the effects of isospin symmetry 
breaking being considered [37H39{ . Consequently, it will make it possible to study the properties of neutron stars and 
the structure of neutron-rich nuclei within the present framework. The corresponding works are under way. 
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